Theory for Superconductivity in Iron Pnictides at Large Coulomb U Limit 
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Superconductivity in iron pnictides is studied by using a two-orbital Hubbard model in the large 
U limit. The Coulomb repulsion induces an orbital-dependent pairing between charge carriers. The 
pairing is found mainly from the scattering within the same Fermi pocket. The inter-pocket pair 
scatterings determine the symmetry of the superconductivity, which is extended s-wave at small 
Hund's coupling, and d-wave at large Hund's coupling and large U. The former is consistent with 
recent experiments of ARPES and Andreev reflection spectroscope. 
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Superconducting (SC) iron pnictides have the highest 
transition temperature next to the cupratesfH 0, d, 
H 0]- The parent compounds are metallic spin density 
wave (SDW) state @, 0, 1(J 11 1 . Superconductivity oc- 
curs when part of Fe 2+ ions are replaced by Fe + . A multi- 
orbital Hubbard model may be a starting point to stud> 



may De a starting point to study 
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the superconductivity. 

20l . 21, Since the parent compound is metallic, most 
theories examine the SC instability from weak Coulomb 
interaction point of view EQEEHEl- On the 
other hand, the observed magnetic moment in the SDW 
phase is large [23| . indicating importance of spin cou- 
plings. The dynamic mean field theory 12j also suggests 
its closeness to a Mott insulator. This calls for an al- 
ternative approach from the viewpoint of large Coulomb 
repulsion U, which will be the purpose of the present 
letter. 

The electronic states of the compound are predomi- 
nantly Fe-3d orbitals near the Fermi surface (FS) 0, 
111 . |24| , which is comprised of two hole pockets centered 
at r = (0, 0) and two electron pockets at X = (tt, 0) 
and Y = (0, 7r), in the unfolded Brillouin zone (BZ), cor- 
responding to 1 Fe atom per unit cell. Note that the 
buckling of As-atoms reduces the BZ to the square en- 
closed by the dashed lines in Fig. 1. The FS structure 
can be reproduced by a 5-orbital model [3]. The bands 
near the FS are mainly d xz and d yz orbitals [3], and the 
FS in the reduced BZ can be reproduced by a 2-orbital 
model, which shifts a hole Fermi pocket from the T- to 
the M = (tt, 7r)-points in the unfolded BZ. In this letter, 
we use the 2-orbital model to study the superconductiv- 
ity at large U limit. We argue that our qualitative results 
will remain unchanged due to the simplification of the 2- 
orbital model. We find that the virtual hopping induces 
orbital dependent pairings of charge carriers. The intra 
Fermi pocket pair scattering is strongest, and the pairing 
symmetry is determined by inter pocket pair scatterings 
and is extended s-wave (s±) for small Hund's coupling 



and d-wave for large Hund's coupling and large U. The 
s±- state was proposed by Mazin ct al. [13j based on 
the analysis of the small Fermi pockets and spin fluctu- 
ations, and was found in weak coupling or small U ap- 
proaches^^, EH • Our result appears consistent with the 
ARPES [2{j ] and Andreev reflection spectroscope [2(| . 

The 2-orbital model reads ^ H = H + Hi, where Hi 
is an on-site Coulomb term, and Hq is a tight-binding 
model on a square lattice of Fe- atoms, 
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where e£ m 



is the hopping matrix in k— space, n = 1 or 
2 denote orbitals d xz (or d yz ). fi is the chemical po- 
tential, a = ± represents the electron or upper (+) 
band and the hole or lower (— ) band, corresponding 
to the diagonalized energy £k±- The band and orbital 
representations are related by a unitary transformation, 
Ckno- = J2 a =± u ™( k )cko,T- Here we follow Ref. [20( and 
parameterize Hq by hopping integrals t™ n between two 
sites i and j = i + f , which is the Fourier transform of 



c„ m (k). We set tf = tf = h, if - if = t 2 , t%l y = t 3 , 
and tfj^y = ±t4 by lattice and orbital symmetry. 

By choosing t\ = — t, t 2 = 1.3t, t% = = — 0.85t, 
the calculated FS with electron density per site w 2.10 
is reproduced in Fig. [TJ which is similar to the first prin- 
ciple calculations [13J, |2fJ for LaFeAsO. The weight con- 
tributed from each orbital at the FS is illustrated in the 
figure. The state on the electron pocket around the X 
(Y) is mainly from d yz (d xz ) orbital. The state on the 
hole pocket around the T consists of d yz and d xz orbitals 
equally if k is along the diagonals, and mainly from d xz 
(or d yz ) orbital if along the x or y axis. 

The on-site interaction 

Hi = ^ ] \Ufli m -\fli m ^ + Jc\ m ^cl m ^CffnlCim,'\] 
i;m— 1,2 

+ [Ul2nilafl t 2al + J ~A\aA2ofii\<JlCi2o\ (2) 
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FIG. 1: (Color online) Fermi surface in the unfolded Brillouin 
zone (BZ) of Ho- The square enclosed by the dashed lines is 
the reduced BZ. Color scheme illustrates weights contributed 
from orbitals d xz and d yz . Arrows indicate inter-pocket pair 
scatterings with wavevectors q ~ (0,7r), (tt, 0) and (tt,tv). 
Numerics (positive value: attractive) are the corresponding 
scattering amplitudes A™™(q) in Eq. (7) in unit of t 2 /U at 
J = 0. Not shown is the intra-pocket scattering j4""(q ~ 
(0,0)) = 20t 2 /U. 



where 



U and U12 are the intra- 



and inter-orbital direct Coulomb repulsions, respectively. 
The terms with J are the exchange interaction. By sym- 
metry, U = U12 + 2J. 27] In the limit, U » t, each 
lattice site is doubly occupied in the parent compound. 
Upon electron doping, some sites will have 3 electrons (or 
1 hole). A single hole at site i may interchange with a 
two- hole state at site j, leading to a metallic phase. The 
effective interaction between two single holes on neigh- 
boring sites (i,j) can be derived by using second order 
perturbation theory, and it is given by 
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where rh refers to the conjugate orbital of m, and the 1st 
and 2nd terms in H2 are the pairing interactions in the 
spin singlet and triplet channels, respectively. The spin 
singlet pair operator b nm (ij) = -^(ci n -\c jml - c inl c jm ^), 
and the spin triplet pair operators T$ z can be written 
similarly. In Eq. (3) and formalism hereafter, we use 
hole notation. The results plotted in all the figures, how- 
ever, will be in the electron convention. Castcllani et al. 
271 ] studied the spin-spin coupling for a 2-fold orbital de- 
generate Hubbard model in the context of V2O3. Our ex- 
pression here is equivalent to theirs, although the pairing 



forms were not explicitly given in their formalism. The 
spin triplet states become important at J/U — > 1/3, or 
J — > U12, which can be seen clearly from the term in B. 
Below we focus on the spin-singlet state with even par- 
ity, which is energetically more favorable for J/U not so 
large. The pairing interaction between carriers derived in 
the large U- limit should be relevant to the intermediate 
coupling region (28j |. 

The effective Hamiltonian is then H e ff = Ho + H2, 
subject to the constraint of no more than 2 holes per site. 
This can formally be represented by a Gutzwiller projec- 
tion operator to project all the unphysical states, similar 
to that in the t-J model [2^]. H e ff may be studied by 
using a renormalized Hamiltonian approach to take into 
account the projection [3(3] by introducing renormaliza- 
tion factors, gt for Hq and 52 for H2, both are doping 
dependent. For a given doping, the effect of the renor- 
malization is to scale all the t's to gtt's, and (U, J) to 
(fft /92)(U, J). Below we will absorb these renormaliza- 
tion factors into the parameters (t's and U) and effec- 
tively set g t = 52 = 1 in our calculations. 
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FIG. 2: (Color online) Energy per site of H e ff in s± state 
(red solid line) and d-wave state (blue dashed line) for (a): 
t/U = 0.1 and (b): t/U = 0.2. 

H e ff can then be solved using a mean field theory 
by introducing mean fields for the spin-singlet pairing 
with even parity and symmetric orbitals [3l| . A nm (r) = 

A m „(f) = ^(b nm (i,i + T)), withf = ±x, ±y, ±(x± 

y). By symmetry, depending on s± (Ai g ) or d-wave (Bi g ) 
states, we have A n (x) = ±A 22 (y), A n (y) = ±A 22 (£), 
Ai 2 (x) = A 12 (y) = 0, An (5 ± y) = ±A 22 {y T x), 
A i2 (i + y) = -A 12 (i - y). Note that A 12 (x ± y) = 
for the d-wave state. The pairing strength with A2 g and 
B2 g symmetries [3l| ar e found very tiny, and will not 
be discussed further. [32] The mean field Hamiltonian of 
H e ff can be written as 
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where = f^+f ^~ k +l' S -k-i)- V(k) is a 2 x 2 
matrix in band picture, given by 

V a0 (k)= E A^\f)A* nm (r)^u m , a (\,)u n ,^) 

nmm'n' ' jr 
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Hmf can be solved self-consistently, and the energy per state, |V++| >> V+_|, \V— \ on Y-pocket. This suggests 



site is E — — ± E± (k) , with E± (k) the quasi- 

particle energy of the upper (+) and lower (— ) bands, 
given by 



E ± (k) 



+ V\_ ± x /wi + V£_ [(SO 2 + 4V 2 ] (6) 



where 5£ = £ + V" = [V++ + Vl_]/2, and u£ = [^ + 
| ± + K 2 _)]/2, and the k-dependence is implied. 
In Fig. [2j the energies of the SC states are depicted as 
functions of J/U for t/U = 0.1 and t/U = 0.2. At t/U = 
0.2, the s± state is always energetically favorable. At 
t/U = 0.1, the ground state is s±-wave if J/U < 0.16 
and a d-wave if J/U > 0.16. 
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that the SC pairing is mainly due to the electron pairing 

of the same orbital. At the pocket centered at T, V-\ 

is negligibly small, so that the SC pairing is mainly due 
to the hole pairings. We emphasize that although there 
are nodal lines, V++ on pocket Y and V— on pocket T 
are always finite. The quasi-particle energy on the Fermi 
pockets are given by E- (k) , which are shown in Fig. |4je) . 
There is a full gap on both Fermi pockets around Y and 
F, consistent with recent ARPES and Andreev reflection 
spectroscope results. Because of the above analyses, we 
have £_(k) w V++(k) around Y and B_(k) w V—{k) 
around T. The results for the d-wave state are also shown 

in Fig.[U The nodal line of V (k) crosses the hole Fermi 

pocket and leads to a d-wave like quasiparticle spectrum 
. The quasiparticle energy at the nodal point is given by 
£ k = V 2 _(k)/E + (k). Since F+_(k) ^ 0, but small, E k 
is non-zero but very tiny [not distinguishable from in 
Fig. 0(f)]. 
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FIG. 3: (Color online) Intra-band pairing amplitude V a , a for 
s± (left) and d-wave (right) symmetry states. Upper panels 
(a,b): electron band, middle panels (c,d): hole band. Fermi 
surfaces are indicated by the black lines. Lower panels (e,f): 
relative sign of the pairing amplitudes and order parameters 
around Fermi pockets. 

In Fig. 02 we plot the intra-band pairing amplitude 
V+ + (k) for the electron-band and y__(k) for the hole 
band. In the s± state, V(k) is invariant under a 7r/2 
rotation, and V++(k) and F__(k) have a nodal line in 
the BZ. V++ have the same sign on X and Y pockets, but 
are opposite to y__ on T. In the d-wave state, V aa (k) 
changes a sign under a 7r/2 rotation, and has nodal lines 
along the diagonals in the BZ. 

Let us examine the pairing strength at the FS around 
the Fermi pockets Y and T. For a Fermi wavevector 
kp on the Fermi pocket centered at C = (fe^fcp, we 
define an angle 9 = arctan(fc^ — ky)/(k£ — The 0- 
dependences of V(k) are plotted in Fig. 0] For the s± 



FIG. 4: (Color online) Angle dependence of pairing amplitude 

V++(k) (red line), V__(k) (blue line), and V-\ (k) (green line) 

along the Fermi pocket around Y (T) in the s±-state [panel 
(a/c)] and the d-wave state [panel (b/d)]. (e) and (f): the 
quasiparticle gap on the electron Fermi pocket (red line) and 
hole Fermi pocket (blue line) for s±- and d-wave states. 

To better understand the SC pairing and its symme- 
try found above, we examine the pair scatterings in the 
orbital representation (intra and inter-orbitals) near the 
Fermi pockets. The spin singlet pairing interaction in Hi 
can be written as 

R Z h = -^ E ^™ m "'(q)&L(k)fo"' m '(k'), (7) 

kk',nmm'n' 

with q = k k', b nm (k) = ^(c k „ T c_ kmi - c klU c_ kmT ). 

4™ m "'(q) is the Fourier transform of A™ r r /(ij). H% rb de- 
scribes the pair scattering processes between two pairs of 
electrons with momentum (k, — k) and (k',— k'). Much 
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of physics may be gained by examining the orbital di- 
agonal term b nn (k). Denote A nn /(q) = (q), with 
q = (qx,q y ), we find 

~ , , 4U(t\c x + t 2 2 c v ) . 1 2 x< 2 

yH ' U 2 -J 2 U + J u-j' 3 y 

Mi (q) = 

^12 (q) = ^ , j l \ c * c v ~ jj2_ ja ^icx + Cy) ( 8 ) 

where c K = cosg^, c y = cosq y , and we have set t± = t 3 
for simplicity. Since we have small Fermi pockets, the 
pair scattering wave vectors are q rj (0, 0) within the 
same pocket , and q « (it, 0) or (0, it) between the pock- 
ets r and X or Y, and q w (7r, 7r) between the pockets 
X and Y, as illustrated in FigfTJ From Eq. (8), we 
find that the intra-pocket pair scatterings are always at- 
tractive (A(0, 0) > 0), and strongest between the same 
orbital, and the pair scatterings between hole and elec- 
tron pockets are always repulsive (A(0,n) < 0). The 
pair scattering between the two electron pockets at X 
and Y points is mainly between two different orbitals, 
and v4i2(7r, 7r) is attractive at small J/U, and repulsive 
at large J/U. This qualitatively explains the relative 
signs in the order parameters among the different Fermi 
pockets in both s± and d-wave states as shown in Fig|3] 
The scattering amplitudes in the case J = are shown 
in Fig. 1, which is of s±-symmetry. 

We have used Eq. ([7]) and © to examine the effect 
to the superconductivity due to the simplification of the 
2-orbital model, which results in the shift of a hole Fermi 
pocket from the T- to M- point. We have found that 
the qualitative physics obtained from our study of the 
2-orbital model remains the same except the parameter 
space for the extended s-wave state is enlarged when more 
accurate band structure is considered. To further ensure 
the qualitative conclusions of our theory, we have exam- 
ined a 3-orbital model as in Ref. [Til], in which there 
are two hole pockets around T in the unfolded BZ, which 
is better in agreement with the LDA calculations. We 
have extended our analyses of Eq. (7) to that model and 
the pairing symmetries are found essentially the same as 
from the 2-orbital model. 

In summary we have examined superconductivity in 
iron pnictides using a 2-orbital Hubbard model at the 
large U limit. An extended s-wave pairing is found most 
stable in a large parameter space, consistent with early 
theories starting with weak coupling (small U) and with 
ARPES [25| and tunneling experiments [2^]. Contrary 
to some of weak coupling theories, we find that the pair- 



ing is mainly from the pair scattering within the same 
Fermi pocket. Our analyses suggest some similarities be- 
tween the superconductivity in iron pnictides and in the 
cuprates. We wish to acknowledge the partial support 
from RGC grant of HKSAR and from Swiss National 
Foundation through the MANEP network. 
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